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We discuss optical absorption in topological insulators and study possible photoelectric effects 
theoretically. We found that absorption of circularly polarized electromagnetic waves in two- 
dimensional topological insulators results in electric current in the conducting ID edge channels, 
the direction of the current being determined by the light polarization. We suggest two ways of 
inducing such a current; due to magnetic dipole electron transitions stimulated by irradiation of 
frequency below the bulk energy gap, and due to electric dipole transitions in the bulk at frequencies 
larger than the energy gap with subsequent capture of the photogenerated carriers on conducting 
edge states. 

PACS numbers: 73.25.+i,72.25.Fe, 72.40.+W 



The rich new physics and unusual properties of topo- 
logical insulators (TI) attract much attention in recent 
years, for review see [ll-Q. The spectrum of TI in the 
bulk has a finite energy gap Eq, while there are topo- 
logically protected conductive surface states inside the 
gap. The states with somewhat similar properties were 
predicted also long ago at the surface of a conventional 
narrow-gap semiconductor 01 and at the interface of het- 
ero junctions made of two insulators with inverted energy 
bands Q without appealing to topological reasons. The 
chiral states studied in recent works on TI [H-Q and the 
chiral states considered in old papers 0, [Hj were derived 
using different theoretical models and the states have dif- 
ferent spin structure. Derivation of the surface state of TI 
was essentially based on quadratic dispersion relation of 
the mass term in a Dirac-type Hamiltonian which is valid 
for materials like Bi2Se3 l2|, H or HgTe/CdTe semicon- 
ductor quantum wells [ll,|Q]7and the wave functions of the 
surface states are formed from the states of bulk bands 
with the same spin directions. Therefore, the electron 
surface states in TI can be characterized by spin oriented 
perpendicular to the momentum. While in Rcfs. 0, [Hj 
the mass dispersion is not needed to derive the surface 
states, and the surface states arc formed from the bulk 
band states of the opposite spin directions. 

Experimental study of transport properties related to 
the surface states is typically complicated by residual 
electrons and holes in the bulk. In this respect it would be 
interesting to study the possibility to induce the current 
in the surface state by absorption of circularly polarized 
light. As the left/right circular polarization of the light 
is associated with the angular momentum projection ±1, 
the absorption of polarized light can vary the number of 
electrons with definite spin. Since the directions of the 
spin and momentum of electrons in the chiral states are 
related, variation of the spin distribution may affect the 
symmetry of the momentum distribution, and, hence, ex- 
cite the current. As these effects are related only to the 
chiral surface state, the bulk contribution to the current 



should be absent. 

Below we set e, h and fc^ to unity, restoring dimen- 
sional units in final expressions when necessary. 

Semiconductors with few energy bands near the Fermi 
energy are usually described by an approach developed 
in seminal paper by Kohn and Luttinger 0] . In the sim- 
plest case when only two bands are taken into account the 
Hamiltonian has a form of the effective Dirac equation 
for the envelope functions. It has a form of the effective 
Dirac equation with the effective velocity of light c* pro- 
portional to the matrix element of the momentum for the 
interband transitions, and the mass term proportional to 
the energy gap Eg — 2M, 



H = c*(t^ X ak,Pk) + T^M, 



(1) 



where Pauli matrices ak operate in the spin state, and r^. 
refer to the conduction and valence band states. Based 
on such an approach, the chiral interface states with 
Weil spectrum were suggested long ago Q at the surface 
of semiconductors, and at the interface of heterostruc- 
tures consisting of two semiconductors with inverted en- 
ergy bands In the former paper the bulk Hamilto- 
nian was supplemented with phenomenological boundary 
conditions at the surface of the semiconductor derived 
from condition that the Hamiltonian is Hcrmitian in re- 
stricted area and obeys time reversal symmetry. This 
phenomenological boundary conditions were shown Q to 
be consistent with the results for hctcrostructurcs and 
yield the surface chiral states that arc similar to the in- 
terface states in heterojunctions. 

Consider first the chiral surface states described by the 
simplest model [41 that describes the chiral surface states 
depicted in Ref. J] as the Tamm states. The boundary 
condition links the spinors related to the conduction and 
valence bands = mo(cr, n)$c, where n is the normal 
to the surface of the semiconductor, or to the boundary 
of the heterocontact. In case of a free surface of a semi- 
conductor Op is a phenomenological constant the value of 
which is determined by the properties of the surface. In 
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inverted heterocontacts the value of oq was shown [8| to 
depend on work functions of semiconductors forming the 
heterocontact, and Cq = 1 for the symmetric case. 

Eigenstates and energy spectrum Eaip) = Eq — avp, 
V = c*2ao/(l + ao), describing the two-dimensional chiral 
surface states at plane z = are characterized 3 by 
momentum p = {px,Py) = p(cos iy9, sin (^) and by the cone 
number a = ±1, and have a form 



yields 
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Here |c) and \v) are periodic Bloch functions of the ex- 
trema states in the conduction and the valence bands, 
and (x, y) are the coordinates along the surface state. 
Decay of the wave functions in the direction perpendic- 
ular to the surface is described by single exponential de- 
pendence with the exponent A ^ v/Eg. 

Three-dimensional Bi2Se3-type TI is described by the 
effective Hamiltonian [§-ll| similar to ([T]) but with mo- 
mentum depending mass term, M M — Bp^. In con- 
trast to the former case, the boundary of the TI with 
vacuum was usually described by model open boundary 
conditions. The wave-functions of the chiral states with 
the spectrum Saip) ~ —avp are 
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(|P1J) + |P2;))e'P'- (e-^^^-e-^^^). 

(3) 

Here I-Plj) and \P2^ ) refer to orbitals related of the 
conductivity and valence bands [l^l- Decay of the 
surface state in z-direction is described by two exponents 
Ai,2 = 25(^^ ± - AMB). 

In both models the exactly backward scattering by 
non-magnetic defects in the chiral surface states is pro- 
hibited, the probability of clastic scattering with varia- 
tion of momentum directions from the angle ipi to ip2 
being proportional to cos^ ^ However, the elastic 

scattering time is not suppressed much because of the 
sizeable scattering for a finite angle. 

An important difference between chiral Tamm ([2|) and 
TI states ^ is that in the former case the surface states 
are formed from the bulk band states of the opposite 
spin directions, while in the latter case spin directions 
are the same. The finite expectation value of the elec- 
tron spin in the Tamm states ([2]) is non-zero only due 
to asymmetry of the contributions from the conduction 
and the valence bands at ao 7^ 1, s = a[{aQ — l)]/[2(aQ -f- 
1)](— sin (/3, cos (/?, 0). While in the TI chiral states case 
^ the spin expectation value is s = sinip, cosip, 0). 

The difference in the structure of the chiral states re- 
sults in a difference of their optical properties. Calcula- 
tion of matrix elements of vertical optical electric dipole 
transitions from the lower (a = 1) to the upper (a = — 1) 
cone of the spectrum for the Tamm and TI chiral states 
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{c\d\v), dTI = 0, 



where (c|d|?;) is the dipole matrix element for the bulk 
interband transitions. Thus the dipole optical transitions 
between the lower and the upper cones are allowed in the 
Tamm states and are forbidden in the TI states. 

As the electric dipole transitions in the TI states are 
forbidden it is worth to discuss the magnetic dipole tran- 
sitions though they are much weaker. Matrix elements 
of the magnetic dipole transitions between the states ([3]) 
depend on direction of the momentum. 



niT/ cx fiBi^Hx cos Lp,iHy simp, Hz), 



(4) 



where H is the magnetic field of the light. If = 
then for an electron moving at an angle cf) to the mag- 
netic field, according to |mr/p cx sin^ 0. Therefore, 
for linearly polarized light propagating in the direction 
perpendicular to the surface of TI the optical absorp- 
tion in the edge states is anisotropic. This must result 
in anisotropic photoconductivity, but the photoconduc- 
tivity is not expected to be large as optical absorption is 
caused by magnetic dipole transitions. 

Consider now the right (left) circularly polarized elec- 
tromagnetic wave propagating along the y axis parallel 
to the surface and characterized by its polarization. In 
this case the probability of the transitions from the state 
in the lower cone a = 1 to the upper one, a = — 1, is pro- 
portional to (hbHo)^ cos* ^ for the left polarized and to 
(iibHq)^ sin^ ^ for the right polarized light. So there is 
a preferential direction of the momentum of electrons ex- 
cited by circularly polarized radiation. The left polarized 
wave preferentially creates electrons moving in positive 
x-direction and holes moving in the opposite direction. 
Hence, optical absorption must induce an electric cur- 
rent in the direction perpendicular to the wave vector of 
the incident electromagnetic wave. However, we expect 
that this effect is difficult to observe since the circulating 
current is greatly suppressed by elastic scattering. There 
are more chances of experimental observation of the effect 
in two-dimensional TI where the finite angle scattering 
is absent and the backward scattering is forbidden. 

Consider now two-dimensional TI state realized in the 
CdTe-HgTe-CdTe quantum wells [l| and described by the 
effective two-dimensional Hamiltonian for El and HI 
sub-bands Idl 



H = e{p)l2x2 

Hp) = v{pxTx 
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(5) 



We will ignore e(jj) and measure energy from the Dirac 
point since this term just shifts the spectrum while topo- 
logical properties and matrix elements do not depend 
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on e{p). The Hamiltonian is expressed in the basis of 
\E1+), \H1+), \E1-), \H1-) where the states El and HI 
are two sets of Kramers partners originating from s-type 
Fg and p-type Fg bands, respectively. Then the hehcal 
edge states are described by energy dispersion = ^vpy, 
where 7 = + for the branch of right-moving electrons and 
7 = — for the left-moving electrons, and the wave func- 
tions are 

^+it) (X J ^ + z\Hl+))e^Pyf{x), (6) 



\El-)-t\Hl-))e^Pyf{x). (7) 



where f{x) = 



^ . As long as the spins of the 



right-movers and left-movers are the opposite the electric 
dipole transitions between the branches are forbidden. 

Consider a quantum well the edge of which forms a 
closed loop of radius R ^ 1/Ai^2- These conditions is 
not difficult to meet as for a HgTe quantum well of 7 
nm width the decay lengths are estimated as 2 and 30 
nm. In this case we can neglect interaction between the 
edge states at the opposite parts of the sample. We 
choose the z axis perpendicular to the plane of well, 
the X axis perpendicular to the edge and y-coordinate 
along the edge and consider magnetic dipole transitions 
induced by the magnetic field of the circularly polarized 
electromagnetic wave of frequency w < 2|M| propagat- 
ing along z-direction. For the right (left) polarization 
r = 1, (—1) the probability of the transitions from the 
left-mover state e_ to the right-mover state e+ is de- 
scribed by the matrix element of the magnetic dipole 
transition m = [LbEIq(\ — r)/2. Thus optical transitions 
are induced only by left polarized waves which excites 
right-moving electrons and left-moving holes and, hence, 
results in current along the edge states. 

The intensity of magnetic dipole transitions is very 
small, but probability of the radiative recombination is 
determined by the same small matrix element so that 
such a recombination can be quite effective. Some pro- 
cesses which are not taken into account in the model of 
TI (say, due residual magnetic impurities or due to tran- 
sitions involving remote energy bands) may cause an ad- 
ditional recombination, but as the initial and final states 
has the opposite spins such transitions are also very weak 
if there are no magnetic impurities. We will take into 
account such processes phcnomcnologically by means of 
recombination time tr. 

Distribution functions n^{p) of right- moving and left- 
moving electrons are determined by coupled kinetic equa- 
tions. We derive the equations using Keldysh diagram 
technique similar to the original derivation 1^ . Interac- 
tion with electromagnetic fields we take into account by 
means of the standard electron-photon Hamiltonian. As 
the kinetic equations have a transparent physical mean- 
ing and can be understood on phenomcnological grounds 



we do not concentrate on a strict derivation. The kinetic 
equations in a uniform case have a form 

dn~. 



1{G-R) 



/e{n^}. 



(8) 



dt tr 

Here G and R describe radiative generation and recombi- 
nation related to magnetic dipole transitions. The term 
with the recombination time tr in (|8]) phenomenologi- 
cally takes into account additional recombination men- 
tioned above, np{e) is the Fermi distribution function, 
le is the inelastic collision integral describing the energy 
relaxation. 

Probability of optical transitions is determined by 
Fermi golden rule, and the generation rate is given by 

G ^2TTfi%H^{n^ ~n+)S{2pv ~Lj), (9) 

were we neglected photon wave vector k p(v/c) com- 
pared to electron wave vectors. 

To derive the recombination term we have related mag- 
netic field to the field operators of the photon field. Then 
after averaging over photon field distribution we obtain 
Green's function of the electromagnetic field and 
treating them as the equilibrium ones we obtain finally 
the recombination term of the form 



R = A 



ck 



(27r)3 ZkijJ 

[(1 + N^)n+{1 - n_) - N^n_{l - n+)], (10) 

where N^^ is the Planck distribution function and k, is an 
ambient dielectric constant. 

We consider the uniform and the steady-state case, so 
that the right-hand side of Eq. ([8]) vanishes. Then we 
assume that the life time of the photoexcited carriers is 
much larger than the energy relaxation time. This is 
very realistic assumption that is often fulfilled in semi- 
conductors at low temperatures. Moreover, in our case 
the lifetime is much larger than in conventional semi- 
conductors because the recombination needs a spin flip, 
while the energy relaxation within the same branch can 
be quite effective due to interaction of ID electrons with 
3D phonons. Such a process is not limited by conserva- 
tion laws as it would be in the system with ID acous- 
tic phonons. Therefore, in the leading approximation 
the solution of the kinetic equation is satisfied by the 
Fermi distribution functions but with a quasi Fermi level, 
= npis — fJ"y)- Conservation of the total electron 
charge yields relation between quasi Fermi levels of elec- 
trons in different branches, /i+ -I- /i_ = 0. Then we in- 
tegrate equation ^ over momentum. The leading term, 
the inelastic collision term vanishes after integration be- 
cause of the conservation of particle number. In the lead- 
ing approximation we can insert = npi^ — fJ"y) in the 
other terms. Thus we obtain the balance equation which 
at low enough temperatures T uj < \M\ has a form 
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(11) 



where the second term in the left-hand side describes the 
radiative recombination, and the magnetic field of the 
light is expressed in terms of the radiation intensity as 

Then according to Eq. (|lip we find that the quasi 
Fermi level of right-movers increases with intensity W , 
and at W « Wq starts to saturate in the region of ener- 
gies of the order T when reaches the value /i+ = lo /2—eF- 
If we do not take into account processes forbidden in the 
TI model and described by tr then in dimensional units 



.3/2 



For typical values of huj — 2eF ~ 10 meV we get Wo ~ 0.1 
W/m^. Additional recombination described by tr would 
increase Wq. 

The excess density of ID carriers equals SN+ = 
—SN- = /x_|_/(27ri'7i), so the current induced by the pho- 
togenerated carriers is 

/ = ev{6N+ - 6N^) = 2GoAi+/e, 

where Go = e'^/h is the conductance quantum. The 
current increases with radiation intensity increasing (see 
Fig. [1]) and saturates a.t W > Wq at the maximum value 



For hio — 2£i? « 10 meV one can estimate the maximum 
current as / w 1 [iK. 




FIG. 1. Qualitative sketch of the edge current as func- 
tion of intensity of the circularly polarized light of frequency 
cj < 2|M| for the case when the radiative recombination dom- 
inates. 

We considered the case of the left-polarized radiation 
resulting in transitions from left-movers to right-movers. 
Applying similar arguments to the right-polarized radia- 
tion we find that it results in transitions from the branch 



of the right-movers to the branch of the left-movers and 
induces the current in the opposite direction. 

The electronic structure of the CdTe-HgTe-CdTe quan- 
tum wells gives an additional opportunity to induce the 
current in the chiral edge states by circular polarized 
light. The lower 2D band E\ is of s-type and contains two 
components of angular momentum projections mj = 1/2 
and TO J = —1/2, while i?l band consists of two p-type 
components with projections mj = 3/2 and raj = —3/2. 
The symmetry of the bands is similar to that of the con- 
duction and the heavy-hole bands in GaAs. Therefore, 
we can ap ply to our case the arguments by Dyakonov 
and Perel IJ] and conclude that the circular polarized 
light of frequency w > 2\M\ with fc- vector in z-direction 
must induce the spin orientation of photocxcitcd elec- 
trons. Namely, the right polarized light leads to electric 
dipole transitions of electrons from mj = 1/2 states of 
the lower band to toj = 3/2 states of the upper band ex- 
citing electrons with spin up. As the excess electrons has 
spin up they can be captured in the ID branch only of the 
right-moving electrons with spin up and such electrons 
will contribute to the current along the edge of the quan- 
tum well. There will be also photoexcited holes compen- 
sating the current. But one can separate photoexcited 
electrons and holes in space by means of the gate voltage 
applied to induce a potential well for the holes in the cen- 
ter of the specimen and the potential well for electrons 
near the edge of the structure, as it is shown schemati- 
cally in Fig. [2] Then electrons can reach the edge states. 





FIG. 2. A model structure to excite the current along the edge 
states by absorption of circular polarized light of frequency 
cj > 2|M|. ID edge spectrum is shown schematically, the 
states filled by electrons are marked by dense lines, vertical 
arrows designate spin directions. 

but holes cannot, and after energy relaxation only elec- 
trons fill the edge state. Again, we assume that the en- 
ergy relaxation time is small compared to lifetime of the 
photoexcited carriers. This condition is easy to meet 
since non-equilibrium electrons and holes in the sample 
are separated. So the electrons induce the spin-polarized 
current / = Go{M — eF)/e along the edge state. Note 
that number of electrons captured to the edge states is 
much smaller than the number of photoexcited electrons 
since the ID edge states are a set of measure zero com- 
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pared to 2D states in the quantum well. The maximum 
number of ID electrons, HMl/v ^ 1. for a perimeter of 
the edge states L of the order of few micrometers. 

Such a non-equilibrium state with excess number of 
electrons in the spin-up branch and equilibrium state in 
the spin-down branch is stable provided there are no scat- 
tering between the branches. Again, though there is a 
topological protection from this scattering it may occur 
due to processes which are not taken into account in the 
model, e.g., due to residual magnetic impurities or tran- 
sitions involving remote energy bands. Such processes 
can be characterized by a time Tm which is much longer 
than the energy relaxation time. Then the circulating 
current may decay during time r^. However, if an AC 
gate voltage with frequency smaller than the energy re- 
laxation rate is applied then the AC current along the 
chiral edge states will appear. 

The current along the edge states can be measured 
by its magnetic field. The effect can be detected also by 
means of a tunnel contact between the TI and the second 
electrode. If we subtract from the standard expression for 
the tunneling current the part related to the equilibrium 
distribution functions then we find the contribution to 
the tunnel current due to non-equilibrium part of the 
distribution of ID electrons in the edge states 
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where N2 is the DOS of the second electrode and V is the 
voltage at the tunnel junction. If N2 depends on energy, 
say, the second electrode is a semiconductor with a mini- 
mum of the conductance band corresponding to the band 
gap of the TI, then varying V one can get information on 
Sns{e) from measurements of the tunneling photocurrent 
which must not be equal to zero even if y = 0. 
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